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Abstract. The non-thermal X-rays from the SN 1006 NE rim present char- 
acteristic scale lengths that are interpreted in the context of diffusion of a 
relativistic electron. The adopted theoretical framework is the mathematical 
diffusion in 3D, ID and ID with drift as well as the Monte Carlo random walk 
in ID with drift. The asymmetric random walk with diffusion from a plane can 
explain the scale widths of 0.04 pc upstream and 0.2 pc downstream in the non 
thermal intensity of X-ray emission in SN 1006. A mathematical image of the 
non thermal X-flux from an supernova remnant as well as profiles function of 
the distance from the center can be simulated. This model provides a reason- 
able description of both the limbs and the central region of SN 1006. A new 
method to deduce the magnetic field in supernova remnant is suggested. 

Key words: (ISM:) supernova remnants ; X-rays: general ; 
1 . Introduction 

In order to explain the non-thermal emission from astrophysical objects and the 
energy distribution of cosmic rays, the theory of particle acceleration was developed 
using various types of acceleration mechanisms, of which some are briefly reviewed 
as follows : 

• acceleration from Interstellar Medium ( in the following ISM) irregularities, 
Fermi (1949) 

• magnetic pumping , Parker (1958) 

• strong shocks , Krymskii (1977) , Axford et al. (1978) , Blandford & Ostriker 
(1978) , Bell (1978a) , Bell (1978b) and the detailed review of Drury (1983) 



• scattering by hydro-magnetic waves , Skilling (1975). 

At the same time, the transport equation for the electron energy was set up , see 
for example Parker (1965) , Kirk (1994) , Jokipii (1987) and Berezinskii et al. 
(1990) . 

Further on, recent developments in the quasi-linear theory of particle acceler- 
ation by super- Alfvenic shock waves (see for example Vainio & Schlickeiser (1998) 
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, Klcpach ct al. (2000) and Schlickciscr (2002) for a general review) produced 
detailed results on the spatial diffusion coefficient but did not cover the possibility 
of simulating the X-contours of the astrophysical object. We briefly review the 
effects on electron transport that are due to the dynamic nature of a scatterer , 
see for example Skilling (1975) , equation (9): 

• convection of electrons with the bulk flow of the plasma 

• adiabatic cooling/heating of electrons in regions of diverging/converging flow 

• stochastic acceleration of electrons 

While the last effect may be considered small under some relatively reasonable 
restrictions, the first and second effects are important ingredients in the stan- 
dard transport equation of electrons. The effect of convection would be to cause 
asymmetry in the modelled X-contours by pressing them closer together on the 
upstream side and further away from each other on the downstream side of the 
source. But these two effects require the residence time of relativistic electrons 
in the emitting layer to be greater than the adiabatic/convection time : if this 
inequality is not verified these two effects are negligible. 

The Monte Carlo modelling of particle acceleration and transport is not a new 
concept and has been extensively used in supernova remnant (in the following 
SNR) , in particular: 

• Calculations of test particle spectra and acceleration times are analysed in 
Ellison et al. (1990) by adopting a first-order Fermi particle acceleration for 
parallel shocks with arbitrary flow velocities and compression ratios r up to 
seven, shock velocities ul up to 0.98c, and injection energies ranging from 
thermal to highly super-thermal. 

• Calculations concerning the ion and electron spectra produced by Fermi 
acceleration in a steady state plane parallel modified shock for Mach numbers 
of 170 and 43 are reported in Ellison & Reynolds (1991) . 

• An extended simple model of nonlinear diffusive shock acceleration was de- 
veloped by Ellison ct al. (2000) and includes the injection and acceleration 
of electrons, the production of photons from bremsstrahlung, synchrotron, 
inverse-Compton, and pion-decay processes. 

• X-ray emission in the supernova remnant G347.3-0.5 is explained by an elec- 
tron population generated by diffusive shock acceleration at the remnant 
forward shock, see Ellison ct al. (2001) . 

In astrophysics we see synchrotron emission from places other than the sites 
of active acceleration. It is therefore a question of electron transport and in this 
paper the following questions are posed: 

• The trajectory of the relativistic electron that produces synchrotron radia- 
tion is helicoidal: on which scales can the trajectory be approximated by a 
rectilinear motion? 

• How can a theory of a particle that diffuses through the random walk be 
developed when the mathematical diffusion is adopted? 
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• The rclativistic electrons alone produce e.m. radiation at a given observ- 
ing frequency but the spatial gradients that characterise the X-sources re- 
quire another treatment. Could a theory that simulates the contours of 
non-thermal X-emission be set up? 

• Can a theory of diffusion from a plane in the presence of convection be set 
up? 

• The X-contours of many SNR (for example Cas A, Tycho and SN 1006) 
present a ring enhancement of the intensity toward the external zones: does 
the observed ring have a theoretical explanation and can it be simulated? 

The trajectories adopted in diffusing away from the center of the box were 
described in Sect. 2.3. 

The mathematical diffusion with and without the convection was set up in 
Sect. 3. 

The physical bases of the random walk for relativistic electrons was introduced 
in Sect. 4. 

The Monte Carlo experiment that simulates the diffusion of rclativistic elec- 
trons from a plane in the presence of a stationary state situation was carried out 
in Sect. 5. 

The formula that allows the theoretical counterpart of the intensity of radiation 
to be built was developed in Sect. 6; this mathematical intensity explains the 
observed X-fiux enhancement in the ring of many SNR. 

The two asymmetries in the observed flux are explained with an energy cascade 
function of the field of radial velocities among the various regions of the SNR, see 
Sect. 6.6. 



2. Synchrotron X-ray emission in astrophysical objects 



When X-ray observations from satellites like Einstein, Chandra and ASCA are 
directed toward typical astrophysical objects like SNR , super-bubbles and extra- 
galactic radio-sources, they detect non thermal emission with intensity l(v) oc v~ a 
(here v is the frequency and a the power law index), sec for example Pacholczyk 
(1970) . The SNR, in particular, presents sharp boundaries in the observed flux 
that can be explained by invoking drift effects. On adopting standard values of 
the magnetic field, the gyro-radius of emitting regions turns out to be w 200 times 
smaller than the typical dimensions where the emission is thought to originate. 
This fact allows us to build a theory of electron transport with a step length 2 
decades smaller than the source dimension. 

These observational effects offer a starting point with which 

• to analyse the random walk performed with steps of a physical dimension 
(Jvg/M-i) ' wnere side represents the dimension of the considered box and 
NDIM the number of mean free paths contained in the box, 

• the computation of the synchrotron losses at the light of the random walk, 

• the setting of the radiative transfer equation, 

• to group all the data on SNR in a unique section. 
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2.1. Synchrotron emission and losses 



Once the relativistic electrons are accelerated (see Appendix A for a brief review 
on the acceleration mechanism), the synchrotron losses due to the magnetic field 
are: 

= -Kad E 2 erg sec -1 , (1) 

in cgs E is expressed in ergs , b ra d — 1-579 x 10~ 3 H 2 and H is expressed in Gauss. 
Another type of loss is the adiabatic one, but in Appendix B we have shown that 
they can be neglected because the residence times of the relativistic electrons are 
small. 

We will see in Sect. 4.2.2 that the synchrotron losses are negligible in the X- 
region and consequently the working hypothesis of the random walk without losses 
is justified. 

2.2. Radiative transfer equation 



The transfer equation in the presence of emission only , see for example Rybicki 
& Lightman (1985) or Hjellming (1988) , is 

^ = -k v CI v +j„C , (2) 
as 

where I v is the specific intensity , s is the line of sight , j v the emission coefficient, 
k v a mass absorption coefficient, C the mass density at position s and the index v 
denotes the interested frequency of emission. The solution of equation |2]) is 

J tf (r„) = ji(l-e-^W) , (3) 

where r„ is the optical depth at frequency v 

dr v = k v C^ds . (4) 

The synchrotron emission values of j v Q and k v C, can be found in Hjellming (1988) 
. We now continue analysing the case of optically thin layer in which r v is very 
small ( or k v very small ) and the density Q is substituted with our concentration 
C(s) of relativistic electrons 

j u { = K e C(s) , (5) 

where K e is a constant function of the energy power law index, the magnetic field 
and the frequency of e.m. emission. The intensity is now 

Iv(s) — K e / C(st)dst optically thin layer . (6) 
J s 

The increase in brightness is proportional to the concentration integrated along 
the line of sight. In Monte Carlo experiments the concentration is memorised on 
the grid M. and the intensity is 

= ^ A s x M(i,j, k) optically thin layer , (7) 
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where As is the spatial interval between the various values and the sum is per- 
formed over the interval of existence of the index k. The theoretical flux density 
is then obtained by integrating the intensity at a given frequency over the solid 
angle of the source. In order to deal with the transition to the optically thick case 
in the Monte Carlo simulation , see Sect. 6. 2, the intensity is given by 

I(i,j) = j^^-exp{-K a J2AsxM{i,j,k))) thin i— ► thick , (8) 

where K a is a constant that represents the absorption. Performing a Taylor ex- 
pansion of the last formula ([5]) , the equation {7J is obtained. 

2.3. The trajectories 



Here the approximation of trajectories moving on a straight line is used and 
the randomisation on the two or six main directions is due to the presence of 
scattering clouds (without re-acceleration) , see Appendix A. The length of the 
step 8 is now compared with the length where the standard theories assume that 
the transport is working. If the magnetic field is strongly turbulent, the mean 
free path of an electron can be identified with its gyro-radius. The gyro-radius 
of the electron can be parametrised by using the synchrotron emission frequency 
expressed in MHz units , vm , and the magnetic field expressed in 1CP 4 Gauss 
units, iJ_4i 

P = 2.6 W-^^^pc . (9) 

This should be considered an approximate relationship since even a single electron 
produces an emission spectrum of finite width. On the insertion of the magnetic 
field of SN 1006, ff_ 4 = 0.15 (see discussion on the data of SN 1006 in Sect. 2.4) the 
value of the gyro-radius now depends exclusively on the chosen emission frequency. 
A typical frequency belonging to the radio-window is 1 GHz which corresponds 
to p = 0.141 10~ 6 pc. The thickness of the emitting shell is 0.52 pc which means 
that NDIM should be 3.6 10 6 once we assume that the length of the step is equal 
to the gyro-radius. 

Another zone of the e.m. spectrum analysed here is X-ray emission, see Sect. 6. 
In this case the frequency observed is 4 10 17 Hz which means that p — 2.8 10 -3 pc 
and NDIM=184. 

The transport of relativistic electrons with the length of the step equal to the 
relativistic electron gyro-radius is called Bohm diffusion Bohm et al. (1949) and 
the diffusion coefficient will be frequency (energy) dependent. The assumption of 
the Bohm diffusion allows to fix a one-to-one correspondence between frequency 
(energy) and length of the step of the random walk. 

24. The data on SN1006 

The diameter of the known remnants spans the range from 3 pc to 60 pc and 
attention is fixed on SN 1006 and its possible diameter of 12.7 pc, see Strom (1988) 
. Information on the thickness of emitting layers is contained in a recent study by 
Bamba et al. (2003) which analyses Chandra observations (i.e., synchrotron X- 
rays) from SN 1006. The observations found that sources of non-thermal radiation 
are likely to be thin sheets with a thickness of about 0.04 pc upstream and 0.2 pc 
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downstream of the surface of maximum emission, which coincide with the locations 
of Balmer-line optical emission , see Ellison et al. (1994) . The values of W u and 
Wd are now introduced in the same way as equation (1) by Bamba et al. (2003) 
(the values of x after which the X-ray intensity is smaller by a factor e in the 
upstream and downstream directions respectively) 

f( x ) = J A Cxp -I^l u P stream (1Q) 
^ ' \ A exp — | r£ ^f- 1 downstream . ^ ' 

Here A is the amplitude of the cut and Xq is the position of maximum emission. 
Observation of the six filaments , see Bamba et al. (2003) , gives the following 
averaged values in the case of non-thermal emission: W u = 0.04 pc and Wd— 0.2 pc. 
These observations allow us to build the intensity of an averaged theoretical X- 
filament, Figure Q] 

For reasons that will be clarified in Sect. 3 we introduce b= SNR d ™ meter = 
6.35 pc , a =b -Wd= 6.15 pc and c =b + W u = 6.39 pc. Another interesting 
parameter is the side of the cubical box situated between the internal absorbing 
sphere of radius a and the external absorbing sphere of radius c ; this small box 
will be useful to simulate the transport of particles from the shock region toward 
the upstream and downstream regions. 

The thickness of the layer , A R , is snr radius accor di ng to McCray (1987) 
and allows us to deduce the theoretical value siSe=0.52 pc. Another useful cubical 
box includes all the SNR and has , for example , sidesNR=18-37 pc ; this will be 
called the great box. The value here adopted in the magnetic field , i?4=0.15, is 
similar to H = 10// Gauss of Bamba et al. (2003) ; Dyer et al. (2004) quote a 
smaller value , H = 4 fiGauss . The data presented here plus the spectral index 
on non thermal X-radiation (see Bamba et al. (2003) ) plus the ratio of flux in 
the ring/center regions , Dyer et al. (2004) , are reported in Table [T] together the 
meaning of the symbols. 

One should remember that the profiles of non-thermal X-ray emission in Bamba 
et al. (2003) are found in another diffusion model. They are extracted from the 
total X-ray flux when the thermal component is not negligible. Suggestions for 
a new diffusion model will be presented in Sect. 3 but up to that section the 
exponential profiles in Bamba et al. (2003) will be the observational reference. 
Roughly speaking we can say that given two decreasing profiles of intensity , 
downstream and upstream in respect to the position of the advancing shock , the 
exponential fits can represent a first approximate description of a more complex 
behaviour of the intensity , see Sect. 3. 

3. Mathematical diffusion 



Once the concentration , C, the energy per electron or nucleon , efe , the 
spatial coordinates, r , the diffusion tensor D(r,ek) are introduced and the losses 
are neglected the general diffusion transport equation for a specie (electrons or 
nucleons) is , see for example equation (3.1) in Berezinskii et al. (1990) , 

¥ = v-(% £fc ))vc . (11) 

In general D(r, e^) is variable both for the standard "thermal" diffusion and for the 
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Figure 1: Profile of intensity of a non-thermal filament in the SN 1006 northeast 
shell in the X-band inserted in the small box. Adapted from Figure 4 in Bamba 
ct al. (2003) . 
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Tabic 1: The data on SN 1006 



symbol 



meaning 



vah 



radius of the internal absorbing sphere 
radius of the shock 

radius of the external absorbing sphere 
theoretical distance between internal 
and external absorbing sphere 
side of a big box 
including the SNR 
magnetic field in 10 -4 Gauss 
shock velocity in km s _1 
spectral index(2.0 — lO.OkeV) 



a 
b 
c 

side 

sidesNR 
H 4 

r 

NT nux clnter Tati ° of X - NT flux at limb 

on X — NT flux at center 



6.15 pc 
6.35 pc 
6.39 pc 
0.52 pc 

18.37pc 

0.15 

2600 
(2.1 - 2.4) 
(5.4-10.8) 



Bohm one. Here we will concentrate on the case of constant energy ( or negligible 
losses ) that in our astrophysical case means a given frequency of observation , 
see formula (|53j) . Therefore D(r, tk) is independent from the spatial coordinates 
and transforms in D that is now constant; of course the given energy/frequency of 
astrophysical interest should always be specified because the mean free path (i.e. 
Bohm diffusion ) changes. The diffusion equation (fTTj) becomes the Fick' s second 
equation , see for example Berg (1993) , and it's expression is different according 
to the chosen dimension. In three dimensions it is 

^=DV 2 C , (12) 

in one dimension it is 

m= D ^ ■ (13) 

In one dimension and in the presence of a drift velocity along the x-direction 
u is 

dC = D (P£ _ ^dC 
dt dx 2 dx 

The hypothesis of the steady state allows us to deduce simple solutions for the 
concentration in 3D , ID and ID with drift and for the spectral index. 

The physical justification that allows to set the concentration equal to zero at 
a given distance from the source is reported in Sect. 4.2.3 and Sect. 6.3. 

The solutions for the diffusion equation can also be obtained with the usual 
methods of eigenfunction expansion (see Gustafson (1980) 
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and Morse & Feshbachm (1953) ) and they are cosine series whose coefficients 
decay exponentially; they can be deduced from formula (25) of Ferraro & Zaninetti 
(2004) upon inserting the dimension d. 

3.1. 3D case 

Figure [2] shows a spherical shell source of radius b between a spherical absorber 
of radius a and a spherical absorber of radius c. 

The concentration raises from at r=a to a maximum value C m at r=b and 
then falls again to at r=c . The solution of equation (fT2"j) in presence of steady 
state is 

C(r)=A+- , (15) 
r 

where A and B are determined by the boundary conditions , 

C ab (r) = C m (l - ~) (l - |) ~* a < r < b , (16) 

and 

C bc (r) = C m - l) (| - l) ~* b<r<c . (17) 

These solutions can be found in Berg (1993) or in Crank 1979 . Thus a comparison 
can be made between the observed and the theoretical spherical solution , see 
Figure [3] 

3.2. ID case 



The ID solution of the concentration is the same as the diffusion through a 
plane sheet. If we have a point source at distance b between a point absorber of 
distance a and a point absorber of distance c, the general solution of equation (|13p 
as a function of the variable x that denotes the distance in presence of steady state 
is 

C{x)=A + Bx . (18) 

The boundary condition gives 

C(x) = C m ^— ^ a < x < b , (19) 

o — a 

and 

C(x) = C-J—^- b<x<c . (20) 
— c 

Figure 2] reports the observed X-profiles as well as the profiles of the ID concen- 
tration as given by the diffusion from a plane. 

3.3. ID case with drift 

Also here a point source at distance b between a point absorber at distance 
a and a point absorber at distance c are considered. The general solution of 
equation lfT4"|) in presence of steady state is 



C{x) = A + Be* 



(21) 
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Figure 2: The spherical source inserted in the great box is represented through 
a dashed line, and the two absorbing boundaries with a full line. The observer 
is situated along the x direction, and three lines of sight are indicated. Adapted 
from Figure 3.1 by Berg (1993) . 



X SNR 



11 



t 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1- r 




j , 1 , , , , i , j 1 i 1 . 1 b , 1 1 ■ , d 1 ■ . ■ ■ 1 ■. 

6.15 6.2 6.25 _ „ 6.3 6.35 6.4 



Figure 3: Observed profile of the non-thermal X filament in SN 1006 (full line ) 
and spherical solution of the 3D concentration (dotted line) inserted in the small 
box. Parameters as in Table 1. 
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t 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 f r 




j , 1 , , , , i , j 1 i 1 . 1 b , 1 , ■ , d 1 ■ . ■ ■ 1 ■. 

6.15 6.2 6.25 _ „ 6.3 6.35 6.4 



Figure 4: Observed profile of the non-thermal X filament in SN 1006 (full line 
) and ID solution of the concentration (dashed line) inserted in the small box. 
Parameters as in Table 1. 
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The advection velocity here considered is the downstream velocity of the shock that 
has direction toward the center of the SNR , conversely the triplet a. b, c, determines 
a positive direction toward the advancing shock in the laboratory frame. In our 
application u and X have opposite directions and therefore u is negative ; the 
solution is 

C(x) = A + Be-B x , (22) 
now the velocity u is a scalar. The boundary conditions gives 
C 

C(x) = __ M _, m _ M _ (e~^ x — e"g°) a < x < b downstream side , (23) 
e d — e D a 

and 

C 

C(x) = — g-; m _„ (e~~^ x ~ e~"^ c ) b < x < c upstream side . (24) 
e d° — e o c 

Figure \E\ reports the observed X-profiles of intensity as well as the asymmetric 
profiles of ID concentration as given by the mathematical diffusion from a plane 
in the presence of drift. 

3.4- Variable spectral index in ID 

The diffusion loss equation for relativistic electrons is usually solved by neglect- 
ing the spatial diffusion. The opposite can be done by neglecting the energy losses 
and exploring how the variations in the absorbing boundaries influence the spec- 
tral index. Up to now the more interesting result on the prediction of the spectral 
index resulting from particle acceleration in shocks , i.e. SNR or Super-bubble , is 
due to Bell (1978a) , Bell (1978b) and Longair (1994) . The predicted differential 
energy spectrum of the high energy electrons is 

C(E)dE oc E- 2 dE . (25) 

But this is the spectral index expected where the electrons are accelerated and 
due to the spatial diffusion the spectral index changes. 

Let us consider two energies E\ and E2 (E% > E\) with two corresponding 
numbers of electrons C\ and C2 {C\ > C2) ■ for example with E%/Ei = 10, 

d/c 2 = 100 . 

In the framework of the ID diffusion, the spectral index T of the energy will 
vary along the distance from the injection point according to the following formula: 

T = Log w ^- . (26) 

In order to allow the spatial index to vary, it is enough to modify the absorbing 
distances of the two populations of electrons , a\ —a. and a 2 > a\. This effect 
is reported in Figure O where variation of the spectral index cover the observed 
range in the X-rcgion , see Table 1. 



4. The physics of the diffusion 
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Figure 5: Observed profile of the intensity of non-thermal X filament in SN 1006 
(full line ) and ID solution of concentration in the presence of drift , = 12.2 -pcT 1 
(dotted line), see formula (fBT))) inserted in the small box. Parameters as in Table 1. 
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Figure 6: Behaviour of the spectral index as a function of the distance from 
injection point b in the downstream region of SN1006 ; a\ =a , 02 = 1.0002 a\ and 
other parameters as in Table 1. 
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The reliability of the stationary state, evaluations on the residence times and 
the influence of synchrotron losses on the diffusion are now analysed in the light of 
Monte Carlo diffusion on a discrete ID lattice. In the following, ND1M indicates 
the number of grid points and NPART the number of particles with different 
patterns. 

Stationary state 

It is important to stress that our simulations involving the ID or 3D random 
walk cover the stationary state. Suppose, for example, that an electron reaches 
the boundaries ( one of the two absorbing boundaries ) after n steps. The visi- 
tation/concentration grid and the energy state are also marked for the previous 
(n-1), (n-2). . . 1 steps. This means that the (n-1) step was due to the electron 
emitted one unit of time later than that characterised by n step, the (n-2) step 
was due to the electron emitted two units of time later than that characterised by 
n step and so on . 

This procedure can be replicated for all the n steps ; this is equivalent to saying 
that = or that we are treating the stationary state. 

We now calculate the time scales that the electrons spend inside the simulation 
box in order to test " The stationary state hypothesis" . The time elapsed from the 
injection in the central point is: 

Uw = N step 8(pc)2/c= N step 8(pc) 6.52 yr , (27) 

where N step is the number of steps before reaching one of the two boundaries and 
the time of crossing the distance S is doubled according to formula (f9"Tj) . The 
number of steps necessary to reach the boundaries depends on the numbers NDIM 
and NPART : its distribution is characterised by a maximum value which is 

extracted from the simulation . The inequality characterising the " The stationary 
state hypothesis" is : 

texp ^ trw ; (28) 

where t exp is defined in . 

The stationary state hypothesis can be checked through a numerical and an- 
alytical argument. These two arguments depend fully on the choice of diffusion , 
basically R 2 cx N ; other choices of D can produce different results. 

Here we have assumed that the acceleration times are smaller than the time 
of losses or expansion or in other words the electrons are accelerated in a smaller 
region in respect to that one in which they diffuse. A detailed discussion of the 
concept of the steady state when acceleration and cooling times are comparable 
can be found in Drury (1983) . 

4.1.1. Numerical estimate 



In every simulation the following inequality should be checked: 

N ma X 2 - 5 X tSNR ( r, Q] 



4.1.2. Analytical estimate 
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A simple analytical estimate for the residence time of a particle diffusing out 
from a box of dimension L is presented. The typical diffusing time r is 



2dD 



(30) 



where R 2 (t) is the mean square displacement (see Gould & Tobochnik (1988) 
, equation (8.38 )) , Z)=||<5 is the diffusion coefficient and d is the dimension , in 
our case 1. On assuming that R = L/2 we obtain 

where the mean free path has been replaced by the length of the step 5. The 
steady state hypothesis is represented by 

5 

T <C t exp = -tsNR , (32) 

which corresponds to 

NDIM < t exp = 5 tsNR^- , (33) 
once j has been replaced by NDIM . Assuming L = ^jj" , the inequality becomes 

NDIM « 60 t SNR -^— . (34) 

tiSNR 

2/5 

The radius of a SNR is known to be Rsnr 1 ^ t s ' NR 

V SNR = l^ , (35) 

and the inequality becomes 

NDIM < 24—^ — , (36) 

VsNR 

and given V C N ~ 100 

NDIM < 2400 . (37) 



Residence time 

The existence of the random walk modifies the canonical formula that allows us 
to compute the distance over which a relativistic electron has damped its energy . 
An electron that loses its energy due to radiation damping has a lifetime i> , 

tv w — w SOO^-^-^ec , (38) 
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where E is the energy in ergs , H the magnetic field in Gauss, and P r the total 
radiated power , see Lang (1999) , formula (1.157). The synchrotron radiative 
lifetime of an electron depends on pitch angle ; thus equation (j3"B")) is to be regarded 
as a typical value of r r for an approximately isotropic , highly relativistic electron 
population. 

The energy is connected to the critical frequency, see Lang (1999) , formula 
(1.154) as , 

v c = 6.266 x 10 18 HE 2 Hz . (39) 

The total power radiated in frequency interval , see formula (1.160) in Lang (1999) 
, is given by P(v) oc F(y/v c ) with 

/>OC 

FC/vc) - / K 5/3 ( V )dv , (40) 

where .K5/3 (rj) is a modified Besscl function. The function F(y/v c ) has a max- 
imum at z^/^ c =0.29 and some authors , see for example Reynolds (1998) , use 
this numerical relationship in order to identify the energy of the electron with 
the observed frequency. Is also interesting to compute the averaged frequency of 
emission Vjv c computed as 

v = Ojpwi^ = 

In the following we will use v/v c = 1 that can be considered an average between 
the two points of view previously expressed. 

Upon inserting in (|38|) the energy as given in (|39p : 

r r = 1.251 x X^^^sec . (42) 

The maximum frequency that can escape from the center of the box and the pa- 
rameters for which the gyro-radius equals the length of the step are now computed. 

4.2.1. Maximum frequency allowed 

If a box with a side L has an electron inserted at the center, the classical time 
of crossing the box t c is : 

t c = £ • (43) 



The solution of the equation 



2c 



? = 1 , (44) 



gives the maximum allowed frequency , 



,-rnax = 5 gg 2 x IQ™ Hz , (45) 

L pc Ii_ i 
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where L pc is the side in pc and f/_4 the magnetic field in 10 4 Gauss units . Upon 
inserting the typical parameters of SN 1006 , reported in Table 1: 

y max = q AQ q x 1q 23 Rz _ ^g) 

For greater frequencies the electron cannot escape from the box when the classical 
treatment is adopted . The time t rw necessary to escape from the box in the 
presence of a random walk is 



Uw ~ 4Sv Tr ■ (47) 



The equation to be solved is now 



r- = i • (48) 



By imposing v tr = fj, where Cf is a parameter greater than one, the maximum 
frequency v ™®w carried by the electron is 

max max / 2 2 

"c,rw V c *( Cf NDIM > ■ ^) 

The relativistic electrons move at velocity near to that of the light but the tra- 
jectory is not rectilinear. The parameter c/ averages the irregularities and fixes 
the transport velocity once the length and the direction of the step are fixed. The 
presence of the random walk produces a decrease in the maximum frequency that 
can be carried by the electron . For example, with v™rw = 4-0 10 17 Hz , and 
assuming c/=2 the maximum value allowed for NDIM is 1273; for bigger values of 
NDIM the X-ray emission cannot be sustained. 

4.2.2. Damping length and gyro-radius 

The length travelled by the electron L rw before being damped is 

L 2 rw = 2dDt rw . (50) 
On replacing t rw with it's physical value we obtain the damping length 

L rw = 110.312( ^ ) 1/2 x ( pc . (51) 

Cf NDIM HV*i>l /4 

The electron gyro-radius , see Lang (1999) , formula (1.153), is 

pn2xl0 9 EH- 1 cm , (52) 

with energy in ergs and H in Gauss . It can be expressed in pc by introducing the 
critical frequency ( in Hz) , ( Bohm diffusion) 

1/2 

p«2.67x 10- 19 -^pc . (53) 
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In order to have transport the following inequality should be satisfied 

P < L rw . (54) 
If equation (|53")) is solved for u, the following is obtained: 

v < 3.197 x 10 27 H( ^ f' 3 Hz . (55) 

c f NDIM' v ' 

With the data from Sect. 4, for example, the inequality is verified: the right hand 
side of (JSU) is 6.062 1 20 Hz and the left hand side is v = 4 10 17 Hz. It is also 
possible to test the hypothesis of random walk without appreciable synchrotron 
losses; expressed as an inequality is 

L r w > Lp C . (56) 

Upon inserting the typical values of the X-region (NDIM=183 and v — 4 10 17 Hz) 
formula (j5"Tj) gives 

0.62 >L pc , (57) 

and as being L pc =0.52 , the inequality is verified. In the radio region (NDIM=2 10 7 
and v = 30 10 6 Hz) the application of formula (|5ip gives 

0.7 > L pc , (58) 

and therefore the inequality is verified. 

4.2.3. Number of collisions on an unbounded lattice 



The random walk in ID on an unbounded lattice is characterised by the length 
of the step 6, which in our case is equalised to the value of the electron's gyro- 
radius p, and by the number of steps after which the random motion is stopped. 
The damping time ,r r , is given by equation (|38p and the time ,r p , necessary to 
travel a distance 5 = p at a velocity vt r — -jh is 



r p = 2.66 x 10- 11 V ^ 2 cL/ sec , (59) 

where v c is expressed in Hz and H in Gauss. The number of collisions made by 
the electron before halving it's energy is 

N = — = 4.699 x 10 22 — !— . (60) 

T p V c Cf 

Table [5] reports the set of parameters after which the profile of the concentration 
in the asymmetric ID random walk is similar to the observed one. The number of 
collisions before to damp the relativistic electron is independent from the magnetic 
field but depends from Cf the factor that fixes the transport velocity and v c , the 
frequency of emission. Is also possible to compute the mean square displacement 
R 2 (t) as given by equation (|30p corresponding to the damping time r r 
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Table 2: The physical data on asymmetric ID random walk 

symbol meaning value 

H-4 magnetic field in 10~ 4 units 1.2 

p electron's gyro — radius 1.24 W~ 4 pc 

v c critical frequency in X — region 4 10 17 Hz 

N number of collisions 58744 



Is interesting to point out that the mean square displacement is independent from 
the chosen energy/frequency ( framework of the Bohm diffusion) and therefore the 
observed profiles should be independent from the chosen type of astronomy. This 
independence from the band of observation can explain Figure 4 of Dyer et al. 
(2004) where there are minimum differences between radio and x profiles in the 
image of SN1006. 

The assumption of the mathematical diffusion that sets the concentration equal 
to zero at a given distance from the plane where the diffusing substance is intro- 
duced , is justified by the physical effect here described. 

5. Montecarlo diffusion 

A numerical simulation is now performed by implementing the asymmetric 
random walk of a charged particle on a ID regular network consisting of a lattice 
with NDIM grid points embedded on a length side(pc); the conversion between 
the unit length and physical dimension of one step is the factor (jv^l^-i) • '^ nc 
analysis is divided into random walks with absorbing boundaries and random walks 
on an infinite or unbounded lattice. 

5.1. Absorbing boundaries 

The rules adopted in implementing the ID asymmetric random walk with in- 
jection in the middle of the grid are : 

1. The first of the NPART electrons is chosen. 

2. The random walk of an electron starts in the middle of the grid. The 
probabilities of having one step are p\ in the negative direction (upstream) 
,pi = 2 — asym x i, and p2 in the positive direction (downstream) , p2 = 
| + asym x |. 

3. When the electron reaches one of the two absorbing points , the motion 
starts another time from (2) with a different diffusing pattern. 

4. The number of visits is recorded on M. , a one-dimensional grid. 

5. The random walk terminates when all the NPART electrons are processed. 
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6. For the sake of normalisation the one-dimensional visitation/concentration 
grid M is divided by NPART. 

These transition probabilities can be considered the ID counterpart of the 2D 
formula (4.5) in Ferraro & Zaninetti (2001) . The introduction of the asymmetry 
in the probability means that it is possible to simulate diffusion and convection. 

There is a systematic change of the average particle position along the x- 
direction: 

(dx) = asym 5 , (62) 
per each time step. If the time step is dt= — and the transport velocity 

v tr = — , (63) 

c f 

where c/ is a parameter greater than one , the convection velocity u is 

asym c 



(64) 



The density distribution is actually determined by the ratio u/D rather than 
u or D separately, so using a single u and D throughout the system only means 
that the ratio is fixed to a constant value. Let consider a strong shock moving in 
the laboratory frame with velocity u s . The post-shock speed is 

u, p — u s ^ , (65) 

and the post-shock gas moves in the same direction as the shock. On introducing 
the index u and d , denoting respectively the upstream and the downstream region, 
the ratio of the Larmor-radii is 

^ = |r = 4 > (66) 

Pd ti u 

, see equation (9.8b) and equation (9.12) in McKee (1987) . 

Let now consider a frame in which the shock is at rest . The downstream 
velocity , u dsi is 

Uds = -^u s - u s = - ~u s . (67) 

The downstream velocity is negative and is 1/4 of the shock velocity; of course if 
we take an x-axis with positive direction toward the center of SNR this velocity is 
positive. 

Thus, using a single value for both u (the advection velocity) and D (the 
diffusion coefficient) may be justified, but for these values the downstream velocity 
^j- = 650 km s _1 and the downstream pd (which corresponds to the magnetic field 
of equipartition and the largest radiative region) for 5 should be used. 

The expected value of asym can now be easily deduced from equation (|6"4l by 
inserting u=^ : 

asym = c f * 0.0021 , (68) 
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and Figure [7] reports the observed X-profile as well as the profile of concentration 
as given by the ID random walk with drift. A formula for the ratio -jj , where u 
represents the drift velocity, can be derived in the X-region 

u u NDIM „ , 

— =2-c f ; — = 1.525c, pc _1 . (69 

D c side 



5.2. Unbounded lattice 



The procedure that describes the random walk is the same as the previous 
subsection except for the fact that the random walk at item (3) stops after N 
steps rather than when one of the two limiting boundaries is reached. Figure [5] 
reports the observed X-profiles as well as the profiles of concentration as given by 
the ID random walk with drift according to the physics explored in Sect. 4.2.3. 



6. Theory of the image 



The overall behaviour of SN1006 in the X-region , see for example Dyer et al. 
(2004) , presents a spherical symmetry with a ring enhancement and a central 
depression. When these effects are visualised through a cut in the flux crossing the 
center of the image, a characteristic "U" profile is obtained. These "U" profiles 
can be explained through cuts of the intensity given by the 3D mathematical 
diffusion while the agreement with the theoretical predictions excludes the presence 
of absorption. 

The asymmetry in the observed flux can be simulated through the map of the 
field of velocity . 

6.1. The mathematical image 



The observed flux of the overall morphology of the SNR is presented in the 
form of "U" profile , see dashed line in Figure [5] The theoretical intensity can 
be found by solving the integral equation ( \E§ along the x direction , see Figure [5] 
where three lines of sight are reported. 

The concentrations to be used are formulae (fT6|) and (fT7|) once r = x 2 + y 2 
is imposed ; these two concentrations are inserted in formula ([3]) which represents 
the transfer equation. The geometry of the phenomena fixes three different zones 
(0 — a, a — b, b — c) in the variable y; the first piece , I 1 (y) , is 



r\/b 2 -y 2 py/ c 2 -y 2 

I\y) = / 1C ah dx + / 2C bc dx 

W a2 -v 2 W b2 -v 2 



= 2 



bCm (~\/a 2 - y 2 + a In U/ 'a 2 - y 2 + aj + \jb 2 - y 2 - a In (^Jb 2 - y 2 4 

b — a 

b Cm ^c In ^ \Jb 2 — y 2 + bj — \]b 2 — y 2 — c In ^ y 1 c 2 — y 2 + cj + \J c 2 ~ y 2 ^j 



c-b 



< y < a . (70) 
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Figure 7: Observed profile of the non-thermal X filament in SN f 006 (full line ) 
and concentration of the ID asymmetric random walk (dotted line), NDIM=183 
,NPART=2000 , c/=6 and asym =0.013. Parameters as in Table 1. 
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Figure 8: Observed profile of the non-thermal X filament in SN 1006 (full line 
) and ID asymmetric random walk (dotted line), NDIM=2980,NPART=2000 , 
c/=2, N= 58744 and asym=0.0043. Physical parameters as in Table 2. 
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Figure 9: Cut of the mathematical intensity I crossing the center (full line ) and 
X-ray data (dashed line) extracted from Dyer et al. (2004) . Parameters as in 
Table 1, great box. 
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The second piece , 1 (y) , is 



r\/b 2 -y 2 r\Jc 2 -y 2 

I H {y) = / 2C ab dx+ / 2C bc dx 

JO J y/b 2 -y 2 

bCm (-aln(y 2 ) - 2 ^Jb 2 - y 2 + 2aln [^b 2 - y 2 + bj^j 

b — a 

bCm (cln (^fb 2 - y 2 + bj - ^Jb 2 - y 2 - cln {^J 'c 2 - y 2 + cj + \j c 2 - y 2 ^ 
~ 2 c~b 

a<y<b .(71) 

The third piece , I ln (y) , is 



r\ c 2 -y 2 

I m (y)= / 2C bc dx 
Jo 

bCm (-cln (y 2 ) + 2 cln (^J cP- - y 2 + cj - 2 c 2 - y 2 ^ 
= c~b 

b<y <c . (72) 

The profile of I made by the three pieces ( [TUJ) , ( EI]) and ( [75]), can then 
be plotted as a function of the distance from the center , see Figure or like a 
contours , see Figure fTUl 

The position of the minimum of / is at y = and the position of the maximum 
is situated in the region a < y < b, and more precisely at: 



■J — (b — 2 a + c) a lab — 2 be + ac) 
b — 2 a + c 

The ratio between the theoretical intensity at maximum , I max , as given by 
formula (|73p and at minimum (jj = 0) is given by 

J ., Xiiuicrnlor 



I{y = 0) Denominator 
where 



Numerator 



, (ac + ab-2cb)a\ , / (ac + ab - 2 cb) a\ , (c + b)(b-a) 

-a In — — — c + aln -- — — — 6-2 \ „ c 

b+c-2a J V b+c-2a J V b+c-2a 



i / <«+»> V'' +»''-■ 

b + c — 2 a 



^"i" ": t ""r'' 4w "::'"r'' - 

b+c— 2a I V b+c—2a 
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Figure 10: Contour map of / particularized to simulate the X-emissivity of 
SN 1006. Parameters as in Table 1, great box. 
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6 + c — 2a / V o+c— 2a 



; (a _ c f (c + 6) 

6+c-2a 



and 

Denominator = 

- 2 ac In (a) + 2 6a In (a) - 2 6a In (6) +2 6c In (6) -2 6c In (c) + 2 ac In (c) . (76) 

In order to make the model more realistic the maximum value of / is normalised in 
order to have the same value as the flux in arbitrary units (i.e. 1080) in Figure 4 
from Dyer et al. (2004) . 

From the simulation reported in Figure [5] we obtain 

Jring_ = 7U _ (77) 
^center 

where I r ing represents the maximum value of I ( on the ring) and I ce nter the value 
at the center. The X-observations with ASCA of SN 1006 , see Figure 4 by Dyer 
et al. (2004) , give 

— 1080/200= 5.4 minimum value 



flux cen ter 

- = 1080/100 = 10.8 maximum val 



fluX r i n g 



fluXcei 

here flux r i ng and flux cen t er are the two fluxes in arbitrary units ( 200 is the 
maximum value of the fluxes in the central region and 100 is an approximated 
average value of the fluxes in the central region ) . 

The mathematical diffusion here adopted predicts a value given '' ' ' ' - — 



flux cc7lt< , r 

7.14, which is comprised in the reasonable interval of observations given in equa- 
tion ([75]). 

It is also interesting to note that the CAP and CENTER region models ( 
subsets of the SRESC model) adopted in Figure 4 by Dyer et al. (2004) predict 
, roughly speaking, a central flux that is 

faring 1Q8 Q /m = 37 (?8) 



flux cen te'i 

This value is three times bigger than that observed or, in other words, over-predicts 
the non-thermal emission at the center of the SNR. 

The effect of insertion of a threshold flux , fluxtr, given by the observational 
techniques , is now analysed. The threshold flux can be parametrised to flux max , 
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the maximum value of flux characterising the ring: a typical image with a hole 
is visible in Figure [TT] when fluxtr — flux max /5. The X-ray emission from SNR 
can also be presented through cuts in particular regions, see for example Figure 4 
in Bamba et al. (2003) , where the x-axis is reported in arc-sec and X-intensity 
in counts. We also particularize our results in those units , see Figure H"2l where is 
possible to visualise the position of the maximum of the intensity that is situated 
between a and b. The ID solutions of the concentration in presence of drift rep- 
resented by equation (j23|) and (j24|) cannot be integrated along the line of sight in 
terms of elementary functions and therefore a numerical integration is performed, 
see dotted line in Figure fT2l 

6.2. The case of absorption 

The effect of absorption is easily implemented by adopting formula §5§ and 
fixing the value of K a for which the results of the thin layer approximation are 
reproduced. The value of K a is then gradually increased in order to see how 
the effect of the absorption modifies the "U" profiles in the cut of the intensity 
previously reported. This numerical experiment is visible in Figure [TS] in which 
the transition from a "U" profile corresponding to the thin-thick case is analysed. 

Due to the fact that we do not know if absorption is present or not, all the cuts 
are normalised in order to have the same maximum value (1080). This transition 
explains the anomalous value of y^" 3 '"" 3 = 2.7 of the CAP and CENTER codes 
as due to a high value adopted in the absorption's coefficient. In Figure [13] the 
third (dashed) line reproduces this anomalous profile. We conclude this paragraph 
underlying the fact that the absorption is not relevant in the X-ray frequencies . 

6.3. An unbounded lattice and derivation of the magnetic field 

The case of an unbounded lattice in the presence of asymmetry and the num- 
ber of collisions fixed by the damping length produces profiles in the upstream 
and downstream region that are reported in Figure [8] These profiles can be pa- 
rameterised through exponentials of the type of those reported in equation (flU|) 
characterised by W™ and W™ . This analytical expression of the concentration 
of relativistic electrons can be the base of the numerical integral along the line of 
sight with the algorithm outlined in equation ([7]). The typical behaviour of the 
intensity as a function of the distance from the center is reported in Figure [TJ] . 

In this case the simulation predicts a value given by y^"*""" 8 = 9.61, which 
is comprised in the interval of observations given in equation (|78p . The magnetic 
field adopted in this paragraph is independent from equipartition arguments and 
therefore the algorithm here outlined can be considered a new way to deduce the 
magnetic field. 

6.4- Comparison with other methods 

In a study on the amplification of the magnetic field in SNR, the observed/theoretical 
analysis on the radial profile of RCW86 in the X-ray band ,see Figure 7 in Volk 
et al. (2005) , gives «4.5. 
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10 20 30 40 

arcsec 

Figure 12: Cut of the mathematical intensity / , full line, that covers a spatial 
zoom of Figure O but with a different scaling. The radius of the expanding SNR , 
b , is marked with an arrow. The experimental data are extracted by the author 
from Figure 4 (panel la) in Bamba et al. (2003) . The dotted line represents the 
asymmetric case with = 12.2 pc^ 1 . 
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Figure 13: Cuts of the projection grid / crossing the center: X-ray data extracted 
from Dyer et al. (2004) (dashed thin line ), K a = 10~ 2 (full line ), K a = 1.8 
(dashed ), K a = 10 (dot-dash-dot-dash), K a — 20 (dotted). Parameters as in 
Figure [TU1 great box. Transition from optically thin to thick layer. 
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-6 -4 „ -2 2 4 6 



Figure 14: Cut of the mathematical intensity I crossing the center (full line ) 
and X-ray data (dashed line ) extracted from Dyer et al. (2004) . Parameters as 
in Table 2, c/=2, but H 4 =1.2 . The simulation returns W™ = 1.17 10" 2 pc and 
W r d w = 3.8 10~ 2 pc. 
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In the appendix A of a paper on X-ray synchrotron emission from SNR , Ballet 
(2006) , is possible to find a formula (analogous to our formula (|74|) ) that gives 
the ratio between the brightness at the center of the sphere and at maximum ; his 
Figure Al gives j^ZL ~ 3 ■ 

6.5. The gamma emission 

Up to now we have focused the attention on the X-ray emission from SNR ; the 
gamma ray emission ( 30Af eV < E < 30Gey) is now analysed from two different 
points of view. 

6.5.1. The gamma emission from relativistic electrons 

Taking the lower limit of the gamma emission , E=30 Mev , the corresponding 
critical frequency is v c = 7.44 10 21 Hz and the relativistic gyro-radius with a 
magnetic field of H_4=1.2 , see Sect. 4.2.3, p = 0.017 pc. Due to the short 
time of synchrotron losses , N=l collisions, the scale widths in the upstream and 
downstream direction are equal, W u = Wd — 0.037 pc. In this case the simulation 
predicts a value given by f lux ™ n n — 7.54. 

6.5.2. The gamma emission from cosmic rays 

The cosmic rays with energy range 0.1 GeV < E < 41Q 5 GeV, see Hillas (2005) 
and Wolfcndale (2003) can produce gamma emission from the interaction with the 
target material. Once the CR energy is expressed in 10 15 eV units ( E i5 ) , the 
magnetic field in 10 -6 Gauss (-H-6) the relativistic ion's gyro-radius is : 

Pz = 1.08-^-pc , (79) 

where Z is the atomic number. The maximum energies that allow to apply the 
random walk of CR in shells of SNR are computed. The magnetic field of equipar- 
tition He = 15 , the minimum width of the X-profile , W u = 0.04 pc and the 
proton , Z = 1, are selected. In order to sustain the random walk the following 
inequality should be verified 

Pz < W u , (80) 

that transformed in CR energy is 

E 15 < 0.55 or E < 5.5 10 5 GeV . (81) 

This upper limit is near to the maximum energy of CR that produces gamma 
radiation , 4 10 5 Gev. The various physical processes that allows to produce 
gamma radiation from CR are now summarised : 

• Acceleration of CR at the shock discontinuity up to E ss 5 10 5 GeV, 

• Diffusion of CR with step's length lower than W u , 

• Absorption at W u and Wd due to a change in the magnetic field from 
#4=0.15 in the shell to #4=0.01 outside the shell, 
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• Production of gamma radiation in a way proportional to the concentration 
ofCR. 

6.6. Asymmetric SNR 

The theory of an asymmetric SNR was developed in Sect. 4.1 of Zaninetti 
(2000) in which an expansion surface as a function of a non-homogeneous ISM was 
computed: in the same paper Figure 8 models SN1006. The diffusing algorithm 
adopted here is the 3D random walk from many injection points ( in the following 



1. The first IP is chosen 

2. The first of the NPART electrons is chosen. 

3. The random walk of an electron starts where the selected IP is situated. The 
electron moves in one of the six possible directions. 

4. After N steps the process restarts from (2) 

5. The number of visits is recorded on Ai 3 , a three-dimensional grid. 

6. The random walk terminates when all the NPART electrons are processed. 

7. The process restarts from (1) selecting another IP 

8. For the sake of normalisation the one-dimensional visitation/concentration 
grid M 3 is divided by NPART. 

The IP are randomly selected in space, and the radius is computed by using 
the method of bilinear interpolation on the four grid points that surround the 
selected latitude and longitude, ( Press et al. (1992) ). The radius will be the 
selected value + i?/24 in order to generate the IP where the action of the shock is 
maximum. 

Our model gives radial velocities , V t heo , 2211 km s _1 < V t heo < 3580 km s _1 
and the map of the expansion velocity is reported in Figure \TE\ from which it is 
possible to visualise the differences in the expansion velocities among the various 
regions as well as the overall elliptical shape. 

Before continuing we should recall that in the presence of discrete time steps 
on a 3D lattice the average square radius ,{R 2 (N)), after N steps (see Gould & 
Tobochnik (1988) , equation (12.5 )) is 



from which the diffusion coefficient , D , is derived 



The two boundaries in which the random walk is taking place are now represented 
by two irregular surfaces. It is possible to simulate them by stopping the random 
walk after a number of iterations N given by 



IP) 



(R 2 (N)} ~ 6DN 



(82) 
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Figure 15: Map of the expansion velocity relative to the simulation of SNR1006 
when 190000 random points are selected on the surface. The physical parameters 
are the same as in Figure 8 of Zaninetti (2000) . 
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where R pc represents the averaged radius in pc . These are the iterations after 
which according, to formula (183|) . the walkers reach the boundaries at a radial 

distance given by from the place of injection; in other words we are working 
on an unbounded lattice . The influence of velocity on the flux F of radiation can be 
inferred from the suspected dependence when non-thermal emission is considered, 
see equation (9.29) in McKee (1987) , 

F = Xt^Hn vl , (85) 

where xt represents the efficiency of conversion of the unitarian flux of kinetic 
energy, fijj the mass of the hydrogen nucleus, no the particles/cm 3 and v s the 
velocity of the shock. 

Assuming that the flux reversed in the non-thermal emission follows a similar 
law through the parameter xx ( the efficiency in the X-region) the effect of velocity 
is simulated through the following algorithm. Once the IP are spatially generated, 
the number of times NTIMES over which to repeat the cycle is given by 

NTIMES = 1 + NTIMESmax * ( v ~ Vmm )3 ; ( 86 ) 

Umax ^min 

where NTIMESmax is the maximum of the allowed values of NTIMES minus 1, 
and v is the velocity associated to each IP. The asymmetric contour map obtained 
when the spatial step is ~ 2 * gyro-radius is reported in Figure \TE\ and the cut 
along two perpendicular lines of the projection grid in Figure ITTl 

In Figurc fTTl thc asymmetry both in the peak to peak distance and the difference 
in the two maximum is evident. The ratio between the X-ray emission in the bright 
limbs (NE or SW) and toward the northwest or southeast ( at 2 keV) is around 10 , 
see Figure 5 top right in Rothenflug et al. (2004) . Conversely our theoretical ratio 
, see Figure [T71 is 9.84. It is also possible to plot the maximum of the theoretical 
intensity as function of the position angle , see Figure 1181 The reader can make a 
comparison with the observational counterpart represented by Figure 5 top right, 
dashed line in Rothenflug et al. (2004) . 

7. Conclusions 

The results of astrophysical interest are now summarised: 

1. A careful analysis of X- intensity of SNR near the external regions predicts a 
complex behaviour as given by the three equations ( 17011 . flTT ]) and ( 172"]) ): the 
fit with exponentials as adopted in Bamba et al. (2003) , our equation (fTQl . 
is not supported by theoretical arguments. 

2. The overall intensity profiles of SNR in the X-region, (represented by the 
new formulae ( [70]) . ( I7T|) and ( 172"]) ) show a characteristic "U" shape that 
can be explained by the mathematical 3D diffusion. The predicted ratio 
between the intensity in the ring region and in the center region (7.15) is in 
agreement with the observations (5.4 -10.3). 

3. The agreement between observed and theoretical overall "U" profiles allow 
us to exclude the presence of self-absorption (optically thick case), see sim- 
ulation in Sect. 6.2. 



X SNR 



39 




40 



L. Zaninetti 




Figure 17: Two cut along perpendicular lines of / reported in FigureflUl. Optically 
thin layer. 
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Figure 18: Azimuthal maximum of the profiles of intensity as function of the 
position angle in degrees. Same parameters as in Figure [TBI 
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4. The Monte Carlo ID random walk with length of the step as given by the 
Bohm diffusion provides the connection with the physical processes as well 
a new formula for the damping length (|51[) . for the mean square displace- 
ment (|61|) and for the diffusion coefficient in the presence of drift (|6"9"f . The 
mean square displacement , equation (|6"Tj) . turns out to be independent from 
the selected band of synchrotron emission . This independence explains the 
morphological similarities between the various regions of the e.m. spectrum 
in SNRs. 

5. The SNR analysed here ,SN1006, presents a spatial asymmetry at the northwest- 
southeast axis as well as a corresponding asymmetry in the flux. They can 
both be simulated once the concepts of diffusion from many injection points 
and variable flux of velocity are introduced, see Sect. 6.6. 

The asymmetric random walk on an unbounded lattice in ID once coupled 
with the concept of damping length can provide : 

• a new way to deduce the magnetic field , H-4—1.2 , against ii_4=0.15 given 
by the equipartition arguments 

• a reasonable explanation for the concentrations of relativistic electrons 

• a value for the ratio between the intensity in the ring region and in the center 
region (6.88) which is in agreement with the observations (5.4 -10.3). 

The analysis in the gamma region can be carried out splitting in two the 
analysis : gamma rays from synchrotron emission and gamma rays production in 
the interaction of CR with target material. 

Eqs. (f6"5)) - (|6T)|) are strictly applicable only to non-relativistic plasma and the 
corresponding relativistic equations based on the first analyses of the theory of 
relativistic Brownian motion, can be considered a target for the future, see Dunkel 
& Hanggi (2005a) , Dunkel & Hanggi (2005b) , Dunkel & Hanggi (2006) and 
Dieckmann et al. (2006) . A subject of further investigation is the role that 
the acceleration of electrons and cosmic rays plays on the shock jump conditions. 
According to McKee (1987) the key parameter is 



w 



Pri 
Ps 



(87) 



where p rs is the post-shock pressure in relativistic particles and p s is the post- 
shock pressure. The effect of the flow of energy from the shock to the relativistic 
particles is to reduce the temperature in the post-shock region and to modify the 
structure of the shock front. 

ACKNOWLEDGEMENTS. The numerical code uses the subroutines extracted 
from "NUMERICAL RECIPES" for mathematical help. The plotting packages are 
PGPLOT 5.2 developed by T.J.Pearson and PGXTAL developed by D.S. Sivia. 



Appendix A. Particle acceleration 

The concept of stochastic acceleration started with the work by Fermi (1949) 
in which it was found that the collisions of a particle of energy E with a cloud 
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produced an average gain in energy (AE) proportional to the second order in u/c 
(u cloud velocity approaching that of light, c ): 

(£>-§(!)■ • m 

This process is named Fermi II. 

On introducing an average length of collision A, the formula becomes: 

dE E 
at t 

where 

It should be remembered that with a mean free path between clouds A the average 
time between collisions with clouds is 

2- , (91) 
c 

see for example Lang (1999) (pag. 467 ) . This means that the particle reaches 
the cloud (our lattice point) in a time that is twice that of straight motion; it is 
therefore possible to speak of pseudo-rectilinear trajectories between the scatterer- 
clouds. 

The probability , P(t), that the particle remains in the reservoir for a period 
greater than t is now introduced, 

P(t)=e-* , (92) 

where T is the time of escape from the considered region. The hypothesis that 
energy is continuously injected in the form of relativistic particles with energy 
Eq at the rate R, produces ( according to Burn (1975)) the following probability 
density ,N, which is a function of the energy: 

, , Rt , E x ^ 
N(E) = —(—)-y , (93) 

with 

7=1 + ^ , (94) 



and r is defined in (|90|) . Equation (|93|) can be written as 

N(E) = K £T 7 , (95) 

where K = E _^y+i ■ A power law spectrum in the particle energy has now been 
obtained; it must be remembered that in this case the ratio of acceleration to 
escape time does not depend on particle energy. 

The strong shock accelerating mechanism was later introduced by Bell (1978a) 
and Bell (1978b) ; the energy gain relative to a particle that is crossing the shock 

■ m 

E 3 c 
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This process is named Fermi I; further on this process produces an energy spectrum 
of electrons of the type: 

N(E)dE oc E- 2 dE . (97) 

This allows us to postulate a kind of universal mechanism that produces syn- 
chrotron radiation with intensity l(v) oc v~ 5 . 

It is interesting to point out that the index 7 of the accelerated particles de- 
pends on the parameters of the shocks u and L , where L represents the thickness 
and u the velocity of the shocks . In particular when u < j- where D is the spa- 
tial diffusion coefficient , the spectrum of accelerated particles becomes very soft, 
see Figure 3 by Ostrowski & Schlickeiser (1996) . The quoted standard spectral 
index of -2 adopted here for the electrons is calculated without escape boundaries 
assuming that there is an infinite box; Ostrowski & Schlickeiser (1996) added the 
possibility of escape and showed its effect on the solution. 



Appendix B. The adiabatic losses 

Other types of losses are those due to adiabatic cooling ,see e.g. Downes et al. 
(2002) , which affects all electron energies and not just the highest ones as in 
synchrotron cooling. Following the approach of Longair (1994) ( equation (19.5) ), 
they are 

~ E^dt>^~^Rdt'~ t exp { ' 

where R is the radius of the SNR and t exp is the characteristic time of expansion. 
The adiabatic losses are negligible if the time scale of the expansion t exp defined 
in equation (|98p is greater than the time , t rw , the relativistic particles have been 
in the emitting shell: 

texp trw • (99) 

In the case of SNR, the analytical solution for the radius can be found in McCray 
(1987) equation (10.27) : 

R(t) = ( 25Ee ^ SNR ) 1/5 , (100) 

where Co is the density of the surrounding medium which is supposed to be con- 
stant , E exp i is the energy of the explosion, and tsNB, is the age of the SNR. 
One should note that this is a Sedov solution, while SN 1006 is young and not 
described by Sedov's formulae. Self-similar solutions for this stage of evolution 
were proposed by Nadezhin (1981), Nadezhin (1985) and Chevalier (1982). Upon 
inserting equation (|100p in equation |9"5|) . the time scale of the expansion takes 
the simple form 

5 

texp = -^tsNR ■ (101) 

Sect. 4.1 shows that with our choice of parameters the inequality given by equa- 
tion ([9"9"|) is always verified; therefore the adiabatic losses can be neglected. 
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